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In this paper, the curvature perturbation generated by the modulated curvaton decay is studied 
by a direct application of <W-formalism. Our method has a sharp contrast with the non-linear 
formalism which may be regarded as an indirect usage of 5iV-formalism. We first show that our 
method can readily reproduce results in previous works of modulation of curvaton. Then we move 
on to calculate the case where the curvaton mass (and hence also the decay rate) is modulated. The 
method can be applied to the calculation of the modulation in the freezeout model, in which the 
heavy species are considered instead of the curvaton. Our method explains curvaton and various 
modulation on an equal footing. 

PACS numbers: 

I. INTRODUCTION 

It is widely believed that a stage of cosmic inflation [l[ which happens in the very early universe is necessary to 
explain the primordial curvature perturbation. During inflation the quantum fluctuations of light scalar fields are 
expanded to become superhorizon classical fluctuations. One or more of those field fluctuations arc supposed to be 
responsible for the curvature perturbation; the field may not be the infiaton field which drives inflation, but the 
curvaton field a 0-111 or the modulation field \ [Mfl, which modulates the decay rate of the infiaton field. 

According to the S N- formalism (see also [l2| ) , the curvature perturbation resulted from a field <f> is given by 

C = SN = N 4 S(f) + ^N^{6<j)) 2 + ^N^(S^) 3 + ■ ■ ■ , (1) 
where N — J din a is the number of c-folds and a denotes the scale factor. The spectrum of the perturbation is 

1 /2 

VA ~ H/2n at the horizon exit. Here 4> may be the conventional infiaton field (ip), the curvaton field (a) or the 
modulation field (x) which are light (compared with the Hubble parameter) during inflation. The fluctuations of 
those fields can (eventually) affect N. The subscript means derivative with respect to <fi. Note that 8<fi is calculated on 
a flat slice (gauge) and C is defined to be the curvature perturbation on uniform energy density slice. It is convenient 
to estimate the magnitude of the higher order effects by using the nonlinear parameters /nl and qml defined by 

C = C 3 + \fNLCl + ^9nlCI + ■■■, (2) 

where Q g = N^Scj) denotes the Gaussian (i.e, the first-order expansion with respect to the Gaussian perturbation 8(f) 
part of (. From Eq. (JXJ) we can see that 1 
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Current experimental data gives (very) roughly [ill [l4{ 

\!nl\ < 100, (5) 
\9nl\ < 10 6 . (6) 



It is also possible to have more than one field (x and a, for example) which can affect N. In this case, besides N atT and N xx we woud 
also have N X<J and the corresponding (6/5)/^ = N X<T /(N x N a ). 
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Po = Po,a + Po, 




a Q (Not Modulated) 



ad (Modulated) 



a c (Not Modulated) 



FIG. 1: The modulated decay in the curvaton model is illustrated. The curvaton density at a is p ,a = \m 2 a 2 <C p a ,-y, which 
gives the curvaton density at the decay pd,a = p ,&X~ :i where X = (a,d/a )- 



In the near future, the PLANCK satellite is expected to reduce the bound to \Jnl\ iS 10 and |<7jvx| ^ 10 5 (HI ^ 
non-gaussianity is not detected 2 . 

In some recent works the question of combiningcurvaton and inhomogencous reheating scenarios has been studied 
in the light of the modulated curvaton decay 19|- 22| , where the non- linear formalism of the component perturbations 
based on , has been used. 

In this paper, we propose a direct method which is conceptually simple and straightforward. In section UH we show 
that the method can be used to produce previous results of [20h22| . in which non- linear formalism has been used. In 
section UlTl in addition to the modulated decay rate, we consider the modulation caused by the modulated curvaton 
mass. In appendix IA 11 we compare the calculation in section |TT] and section IIIII In appendix IA 21 we reproduce 
standard formulas for the conventional curvaton. In section IIV1 the modulation in the freezeout model 0] is solved 
when the massive species might not dominate the density. Our method explains the curvaton model and the various 
modulation scenarios on an equal footing. 



II. MODULATED DECAY RATE OF THE CURVATON MODEL : UNMODULATED MASS 

In this section, we will calculate the curvature perturbation when the curvaton decay is modulated. We consider 
the case in which the decay rate is modulated by a light scalar field \ through a coupling. The case of the modulated 
decay rate through a mass is considered in the next section. 

Consider a curvaton field a with a quadratic potential V(a) = (l/2)m 2 <r 2 . The number of e-folds between curvaton 
oscillation a to a uniform energy density time slice after curvaton decay at a Cl is given by 

N = ln(^) =ln(-) +]n(^)=N 1 +N 2 , (7) 
\a J \a Q J \a d J 

where ad represents curvaton decay at H ~ T and the two parts are denoted by N± and N 2 , respectively. After 
curvaton decay the universe is dominated by radiation therefore H cx a~~ 2 and we can write N 2 as 

where we have used H(a d ) = T and H(a c ) corresponds to the uniform energy density slice (which does not depend 
on a or %). See also FigQ] 

By definition of the oscillating curvaton, when the curvaton start to oscillate some time after inflation when H ~ m 
at a G , the energy density is dominated by radiation 7. The energy densities are hence given by p ~ p D . 7 = 3m 2 M 2 
and p . a = (1/2)to 2 ct 2 . Here the reduced Planck mass is M p = 2.4 x 10 18 GeV and a denotes the field value of the 
curvaton when it starts to oscillate 3 . On the other hand, at curvaton decay (H ~ T) the energy densities are given by 
Pd = 3r 2 Af 2 and pd, a = (l/2)m 2 cr(t ) 2 (a £ i/a ) _3 because the curvaton behaves like cold matter during oscillation. 



2 Note added: PLANCK satellite has recently released their data [l6l. |17H . This has interesting implications for curvaton model in general. 
A curvaton with modulated decay width or mass may help to relax constraints imposed on curvaton scenario by PLANCK data (see 
[lgj as an example). 

3 This is often called g(cr*) where a* denotes the field value of a at horizon exit. In this paper, we only consider the curvaton with a 
quadratic potential, therefore g'(cr*) = dg/da* is a constant which is close to one when slow-roll condition is satisfied. In addition, one 
more derivative makes g"(<j r ) = 0. Therefore we only consider a ~ 0*. The including of g and its derivatives only makes our formulas 
look unnecessarily complicated withour gaining much. 
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If we focus on radiation (which dilutes as a 4 ) during this period and define X = (ad/0,0) for notational simplicity, 
we can write 

3m2 ^ 22 , (9) 
p drt 3r2M2 - sgg' 



3r 2 M 2 X 4 - -m 2 cr 2 X = 3m 2 M 2 . (10) 



For later usage, we define a quantity r used widely in literatures of curvaton by 

ip d ,a 3 (\m 2 o 2 ) X~ 3 S(\m 2 o 2 ) 



ipdn + 3p d ,a 4 [3r 2 M p 2 - ( ±m 2 <7 2 ) X~ 3 ] + 3 (±m 2 a 2 ) X~ 3 12r 2 M 2 X 3 - \m 2 a 2 



(11) 



In the oscillating curvaton model, r is comparable to the ratio of the curvaton energy density to the total energy 
density of the Universe at the curvaton decay. 

The modulated curvaton decay is introduced by the function T(x), which becomes inhomogeneous in space due 
to the modulation caused by a light field \. From Eq. ([T]), the curvature perturbation to linear order is given by 
Q = N x Sx + where other sources (e.g, inflaton perturbation Sip and the conventional curvaton perturbation Sa) 
are included in The last expression in Eq. (jlip is convenient because it allows us to calculate the derivatives of 

r with respect to the fields. A very simple but handy relation which will be used frequently in this paper is 



12r 2 M 2 X 3 _ 12r 2 M 2 X 3 - im 2 a 2 + \m 2 a 2 



r 



12T 2 M 2 X 3 - \m 2 o 2 12r 2 M 2 X 3 - \m 2 o 2 3' 

By making derivative of both sides of Eq. (|TU)) with respect to %■> we obtain 



(12) 



6rr x Af 2 X 4 + 12T 2 M 2 X 3 X X - -m 2 a 2 X x = 0, (13) 



which implies 



6T 2 M 2 X 3 T x 
\m 2 a 2 - 12T 2 M 2 X 3 T 

2 6 r ' 



(14) 



where Eq. (TT2")) has been used. 

N 2x is evaluated by making a derivate of Eq. 



Therefore we have 



N x = N lx + N 2x = ~r 1 f. (16) 

This gives for example, Eq. (31) in [22| and also consistent with [2(| and plj . 

It is also straightforward to calculate the non-linear parameters. By making derivative of Eq. (|16j) once more we 
have 
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In order to evaluate r x = dr/dx, we make derivative of r by using the last expression in Eq. (|llj) and keep in mind 
that r(x) and X{\) are now functions of \. We have 



-3 (§m 2 cr 2 ) 24rr x X 3 + Z6T 2 X 2 X X 
12r 2 X 3 - \m 2 a 2 12r 2 X 3 - \m 2 a 2 



r\ r 



2 11 ^13 11 K| 



3/ r 



r\ X^ 



3/ X 



(18) 



where we have used Eq. (|T2|) to obtain the second equality. From Eq. (fTB]) and (fT4|) we find useful relations 



r x 6 A7 
-f? = —N- 
r 

3 



Y 



X 



x 
X 



+ 1 N- 



(19) 



By using these relations we obtain 



3 - 2r - r 2 



(20) 



Therefore we find 



6 f N xx 
5 fNL 



N 2 r 



- 3-2- 



rr 



XX 



r? 



2-r 



- - 2 - r 



6 /rr 



\X 



(21) 



By taking derivative of Eq. (|A12[) with respect to Xi we have 



2r x 

AT — ± 

O (T 



r(l-r)(3 + r)r x 

9rcr 



(22) 



where Eq. (pO)l has been used to evaluate r x . Therefore by using Eqs. (TIB) and (|A12[) we can obtain 



6 f a X _Nx^_ = (l-r)(3 + r) 



(23) 



We can carry on to calculate gNL- Firstly we substitute Eq. (|2"0")) into Eq. (|T7| to write 



3 - 2r - r 2 Ar2 1 
N 2 - -r 



r xx r - r 2 
r 2 



(24) 



and making derivative with respect to x once more. Then eliminate r x again by Eq. (|20p . finally we obtain 



54 
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135 - 54r - 22r 2 + 10r 3 + 3r 4 - 18 (9 - 2r - 



r 2 

X 



36 



r r 2 

1 XXX L 

r 3 

X 



(25) 



This gives Eq. (43) in [21| and Eq. (28) in [2Cj. It should be clear that our method is capable to reproduce previous 
results. 



III. MODULATED DECAY RATE OF THE CURVATON MODEL : MODULATED MASS 

In this section, we apply our method to the case where the modulated decay rate is sourced by the modulated 
curvaton mass. This is more complicated than the previous case because now the time slice when the curvaton start 
to oscillate is also modulated. In that way, the quantities defined at a D may also depend on x- We thus need to 
introduce a, before a a , since not only but also a Q is modulated. See also FigJ2] 

Before the curvaton oscillation, the curvaton is slow-rolling and we assume that the curvaton density is scaling like 
p a oc a -6 ™ just before the oscillation. In general e w is time-dependent in the radiation-dominated Universe; however 



Pi = Pi,a + Pi 




a, (Not Modulated) 
a a (Modulated) 

a ( i (Modulated) 

a c (Not Modulated) 



FIG. 2: The modulated decay in the curvaton model is illustrated when the curvaton mass is modulated. 

for the modulation at ao we only need e w defined very close to a . Therefore in our calculation we can assume that 
e w is approximately a constant. 

The e-folding number is now given by 



iV = ln(^)=ln(^Uln(^Uln(^ 



N + N!+ N 2 



where N 2 is again given by Eq. © and N 2x is given by Eq. ([15]) , 

In order to calculate No x , we set the fraction of the energy density at ai as 

p Pi,cr _ Pi,<J 



If we define 



we find 



7V = ln(^) ^lnX , 



Pin = 3H?(l-fj) 
3m 2 (l-/ )' 



where 



So = 



r . Fl,o 
Ji — 



Po,a _ Pi, (jXq 

3m 2 M£ ~ 3to 2 M 2 

Pi, a Pi,<r 



Oil - Af 

n — z Ji- 

ox to 

Here, for the curvaton density we assume pi >a oc to 2 . We can find from Eq, (l29[) : 



Considering 



Pi,~r 

Pa,-, 



Xq, we find 



AX 3 X — Pi 2 mxPin I d/oPv-yPo 

° ° X dx Po, 7 m Po.j dx pl n 



^Xqx = 2 m x(P*fi) 2 m x c x Qx (Pofo) 
Xq m p in TO Xq p , 7 



(26) 



(27) 



(28) 



(29) 



(30) 



(31) 



(32) 



4 This is in sharp contrast to the "freezeout and mass-domination model" in Q, which uses p<j oc m 4 for the initial condition. We will 
compare those models in the next section to show that our method explains those different scenarios on an equal footing. 



Solving that equation for , we obtain 



X 2 m \ p in ; 



where 



^Po,*y ~t~ ^wPo,t 



Finally, we have 



(l-O = r-^ ■ (34) 



N 0x = (lnX ) 



l n „ \ m X {-> , P 



For the usual set-ups /j ~ /„ ~ 0, the result shows 



2 m 

For the second stage (i.e, during the curvaton oscillation), we define 



tfox*"^- (36) 



iV"i = In I — ] = In Xi . (37) 



Then 



where 



v4 _ P°,7 _ 3m 2 (1 - f Q ) 

A i - , r9<1 FT) l* 58 ) 

p rf , 7 3r 2 (l-/ d ) 



Pd.cr 


Pi,<jX 1 3 X Q 


Pd 


3r 2 M 2 


% _ 








Po,<7 


Pi.O-Xq 


Po 


3to 2 M 2 


3/o _ 









We thus find from Eq.flM 



, ^ + 2^-2^ I u 

Jo — — 

X /«,. (39) 



4X^ lx = ( 2^ - 2§) ^1 - + d J± P -^2±. (40) 



in 



7 



By using = Xf, we find 



4^ = - 2^ - e ^^M + (_ 3 ** - e w ^x + 2 !!!x _ 2 £x\ . (41) 

Xi \ m T J X a p .~t \ X 1 X m T J p d , 7 



Although not mandatory, we are going to assume e w <C 1 in this paper. Solving the equation for -S^S we find 



X\ x 2(pd,y + Pd,a) /T 



X ""X 



Xi 4pd tl + 3pd,a \ r m 

'1 i\ fv. 



= — -r -+- - , 

6 2/1 r 



\ 



(42) 
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where r is given by Eq. and we have used the relation 



2p d 



4/0 d;7 + 3p d ,cr 4p dn + 3pd,, 



(43) 



Therefore 



N 



(ln*i), 



1 1 
-r + - 
6 2 



Our final result is 



N x = N 0x + N lx + N 2x 
6 r m 



(44) 

(45) 
(46) 



To understand the result, imagine an explicit form of T. Just for instance, one may assume T oc Am' 1 . This gives 5 

1 , , 



-r(l-n)^. 
o m 



(48) 



It might be interesting to note here that when n = 1, the effects of modulated curvaton oscillation and decay cancel 
each other out. 

For the non-linear parameter, we make derivative with respect to x once more to obtain 



N 



1 



1 



vy - -rJl - ti)— + -r(l - n) 
xx 6 XK ' m 6 v ; 



m xx m — m x 



(49) 



where by using similar method as in the previous section, we can obtain 

3 - 2r - r 2 



•_x 

r 



(50) 



Interestingly, r x has the same form as what was given in Eq. (|20[) 6 . The nonlinear parameters are hence given by 

6 



6 „ 3 6 6 m vv m 

t!nl = -r - 2+ n r + —. 7^^- 

5 r r(l — n) r(l — n) m x 



(51) 



and 



54 
25 



N. 



9NL 



x \x 



N 3 

x 



= 3r z + lOr - 4 
18 36 



18 18 

r 



36 
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54 



72 



l — Ti r(l — n) r 2 r 2 (l — n) r 2 (l — n) 2 



54 



108 



1 — 71 r(l — n) 7* 2 (1 — n) r 2 (l—n)' 



36- 



(52) 



IV. FREEZEOUT MODEL WITH THE MODULATED MASS 

The modulated decay scenario of the freezeout model is explained by a massive particle species ip, whose mass M 
and the decay rate T may depend on \. Initially the massive species ij) are subdominant. 



Even if e w is not negligible, we can define n' = (n + e w ) to obtain a simple formula 



1 



(i-«0 



(47) 



This implies also has the same form as Eq. 1 1231) , However they are not really the same because the corresponding N x are different. 
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• In the original model Q it has been assumed that ip were thermalized at some early time (T > Tf > T n ) and 
becomes non-relativistic at the temperature T n ~ M Q . Here Tf is the freezeout temperature. Then the density 
of the massive species ip at that moment is 

Pip(T„) = Mn^,(T n ) ~ M 4 , (53) 

where M denotes the mass of ip species. We usually have the total density p(T n ) > p^,(T n ). Replacement from 
the modulated curvaton scenario is: 

Po,a — 7} m2<j2 Pn,i> K (54) 

which is defined at the beginning of the scaling p^ oc a -3 . 

• Normally, one will assume that the freezeout occurs after ip becomes non-relativistic (M > Tf). In that case we 
find the Boltzmann suppression: 

. . fMT f \ 3/2 ( M\ , , 

MTf) = 9 [^J cxp (-—J , (55) 

where g depends on the model. In that case we find 



3/2 / M 



exp ( -— ) . (56) 



It is possible to introduce a factor [23|, ]24 



TeV (TeV)" 2 ' 
where (av) is the thermal-averaged annihilation cross section, and write 



x f = ^ (57) 

.25 + m^ + m^, (58) 



Pip{Tf) = M A 



9 



(59) 



.(2^/) 3/2 

When the cross section does not depend on M, the replacement from the modulated curvaton scenario becomes 

1 



p 0t „ ~ -m-o- -> pf.^p oc M 3 {x f y 3/2 , (60) 



2- J 



where the trivial identity e ln = M 1 has been used. The modulation about the parameter Xf is weak and 
negligible. The ratio of the energy density at the freezeout temperature is 

(61) 

J1 p f M y 1 

Alternatively, one may assume (av) oc M~ 2 and find 

p M ocM 5 (x f )- 3 / 2 , (62) 

which gives the ratio 

f/E^cc (x f f' 2 M. (63) 
Pf 

Note that Xf = 1 reproduces the first scenario (p^, oc M 4 and / oc M°). In this paper we are assuming instant 
trandition between phases and the simple M-dependence for simplicity. The actual calculation has to be highly 
model-dependent. 
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In this section we mainly consider the first scenario (or cquivalently xt = 1 in the second scenario), since the 
original paper about the freezeout model Q considers that possibility. For the freezeout model, we use the subscript 
"F" to denote the time when the density starts to scale like p^ oc a -3 . If we define 

No = m (2L ) = lnX , (64) 



CI; 



we find 



4 = Pi^ = ft, 7 (65) 

PF. n pf{1 - Jf) 



where pi is defined just before pp so that the density scaling is well approximated by p 7 oc a 4 . We thus find for 

dpi,y/dx ^ 0: 



X 0x 1 dp F 1 df F 



X Pf d\ 1 - If d\ ' 



Therefore 



After the freezeout. we define 



Then, we find 



We thus find 



,x lx _ 1 d PF 1 df F T x , 1 df d 



X 1 p F d x i - If d x r (1 - /„) % ' 



where 



1 dpi? 


1 




Wf 


Pd,ip 


Pf,i/>X 1 3 


Pd 


3r 2 A/2 


dfd_ 




dx 





We find 7 



Solving the equation for -yjS we find 



(66) 



( ^ ) 



^ ~ . (69) 



-2-^ + 7^ FTlP' (70) 



(71) 



2^ ) fd- (72) 



X\ \pf ox r / 1- f F ox 

-^ + —^f t -^) 7T A r Y (74) 



*ix . 2(p d 



Xi Ap d ^ + 3p d 4, T 3pF4, dx \^PF dx J 4(1 - f F ) dx 

6 r + 3p^ 9 X [ '\4pFdxJ i(l-f F )d X 1 J 



' In the simplest case one may assume fp oc M k to find dfp/dx = The first scenario gives k = 0, while the second scenario 

suggests k 0. 
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Therefore 

N lx = QnX 1 ) x 
Our final result is 



r dpF.4, 



(1 



1 dpi 



l-r dfp 



N 0x + N lx + N 2x 



IE* 
6 r 



dpF,4,/dx dp F /dx 



3pF,ip 



4pj 



4p F d X J 4(1 - f F ) 8 X 



ldfp 
4 d x 



(77) 

(78) 
(79) 



To understand the result, consider explicit forms of pF^- 
• For PF,ijj M 4 , k = and pp oc M 4 , we find 

1 



G 



x 

r 



2^ 
M 



(80) 



where r = 1 reproduces the original calculation of ML The result is consistent with the conventional calculation 
of the mass-domination and the freezeout scenario |7[ . 



• Replacing pp^ — > p ocr oc m 2 and pp 
previous section: 



p oc m with k = 0, 8 it reproduces the modulated curvaton in the 

(81) 



1 



N y = r 

x 6 



x 

r 



""X 

rri 



where the e^-dependence does not appear in the above calculation. 

• In a more realistic calculation one must consider Eq. (|56|) and the cross section using the numerical methods, 
which may shift the coefficients [25!]. 



In the multi-field modulation model we find 8N = JA N Xi 6xi + Yli N XiXj 5xi$X 



i r Y 
6 r 
i 



1 dp 



lnT , 
6 3 



3 pF,ip dxi 
\npF,i> - 



r 1 dpp r dfp 
~^PF dxi 4 dx 



1 k 
-\np F + -ln/> 



V-N-y. + r 



— lnT H — lnpp,/, 
6 3 P 



1 /fc 
-lnp F + -In f F 



XiXj 



(r-l)(r + 3) 



N Xi N Xj +r 



ilnT + - ]np F ^ 
6 3 



1 



In pp + - In / F 



(82) 
(83) 
(84) 
(85) 



XiX.i 



In our direct method, it is very easy to evaluate the higher derivatives. 

Notable application of the above result is the conventional curvaton (i.e, the curvaton without extra modulation). 
For the curvaton one needs just a simple replacement Xj ~> a '■ The curvaton hypothesis gives ^f- = and j^;^- — 0. 
Then one will find 



r 1 dp F ,c 
3 pF,a da 



2r 5a 

TIT' 



(86) 



In the above formalism the curvaton mechanism is calculated as a specific example of the modulation. In that way, the 
mixed perturbations of the curvaton and the modulation are calculated in our formalism as the multi-field modulation. 



We find f F ~ p F ,a/pF ~ 2 o- 2 )/(m 2 Af|) oc m°. 
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V. 



CONCLUSION AND DISCUSSION 



In this paper, we proposed a direct application of SN formalism that can be used to calculate the curvature 
perturbation from the curvaton (or the heavy species) with modulation. We calculated for the first time the case 
where the modulated curvaton decay is due to the modulated curvaton mass and obtained non-linear parameters. Our 
method can be compared with the calculation based on the non-linear formalism of the component perturbations. 

Although we consider a quadratic potential for the curvaton which dilutes like cold matter when oscillate, in principle 
the method can be extended to more general cases once we specify the modulation and the dilution behavior [22I l26j ] . 
Our method explains curvaton and various modulation models on an equal footing and provides a convenient way to 
calculate the cosmological perturbations in the multi-component Universe. 

One of the natural cosmological expectations would be that a non-relativistic matter is created and its density starts 
dominating late after reheating. The original curvaton mechanism is based on that simple expectation; however the 
original curvaton mechanism requires significant isocurvature perturbation of the matter density. As a consequence, 
the " non-relativistic matter" is usually replaced by a " sinusoidal oscillation" whose amplitude must be inhomogeneous 
in space. 

In the light of the cosmological model building, the original curvaton conjecture seems to be quite restrictive. There 
could be a deviation from the sinusoidal oscillation (i.e, the scaling of the density could be different due to a deviation 
from the quadratic potential) or the "non-relativistic matter density" could be the "conventional particle" that is 
created by the usual thermal process. 

The former possibility (deviation from the matter scaling) has been discussed in [22j , and the latter (non-relativistic 
"particle" from the conventional thermal process) has been discussed in this paper. The method provided in this paper 
makes the calculation in [22J drastically easy. Obviously, the models discussed in this paper (and in [Hj]) are expanding 
to a great extent the application of the original curvaton mechanism. 

Recently a significant extension of the curvaton scenario has been discussed in [28| , in which an inflationary stage 
is considered for the curvaton mechanism instead of the oscillation. In the name of the "curvaton" , the curvaton 
inflation converts isocurvature perturbations that already exists at the beginning of the secondary inflation into 
curvature perturbations whose wavel eng th is far beyond the reach of the secondary (curvaton) inflation. 9 Higher 
order perturbations are calculated in [22j , although the calculation depends on the indirect method of the non- linear 
formalism. The direct calculation of the SN formalism presented in this paper may have the potential application to 
the inflating curvaton mechanism, which can include any kind of modulation at the same time. 
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We show that by using the same method, we can have a slightly different way to obtain results given in section [TTJ 
This calculation is closer to the one used in section HTT1 

At curvaton oscillating a Q , we set the fraction of the curvaton energy density to be 
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Appendix A: More about the calculation details 



1. Comparison between calculations of Section [TTI and Section IIIII 



Po,a 



(Al) 



3m 2 Af„ 2 ' 



If we define 




(A2) 



9 



Applications of the inflating curvaton can be found in [29H31I ]. In contrast to the conventional curvaton, the non-Gaussianity parameter 
/jvl is expected to be positive in the inflating curvaton, which helps PBH generation in the curvaton mechanism l3ll . 



12 



wc find 



where 



_ Pon _ 3m 2 M 2 (l - f ) 



X i run /»„n 

Pdn 3T*M»(l-f d )' [6) 



U = (A4) 



Therefore 



where p dtl = 3m 2 A/ 2 (l - f )X- i . We thus find 



(A7) 



Finally, we find 

'i i\ r x 



r + -)-£■ (A8) 



tf lx = (ln*) x = -l-r + -j-^. (A9) 
The remaining calculations are the same as those in section [XXJ 

2. Conventional (oscillating) curvaton 

In this appendix, we rederive some familiar formulas of oscillating curvaton [27j by using our method. Let us start 
from Eq. (fTUl) . 

3r 2 A/ 2 A 4 - -m 2 a 2 X = 3m 2 A/2. (A10) 
We make derivative to both sides with respect to a to obtain 

12r 2 A/ 2 A 3 A' - m 2 M 2 aX - \^n 2 a 2 X' = 0, (All) 

which immediately gives 

_ X' _ 2 3(im 2 a 2 A" 3 ) 12 1 

~X ~ 3l2r 2 -imV 2 *~ 3 ^ " 3 r a' (A12) 



The curvature perturbation is given by 



C = N a Sa=lr—, (A13) 
3 er 

which is a standard result of curvaton. 

In order to calculated Jnl, we simply have to make derivative of Eq. (|A12[) with respect with a once more and 
obtain 

2 1 2 1 

N mr = --r- 5 + -r„-. (A14) 

3 <7 Z 3 a 



13 



In order to calculate r a , we have to make derivative of r with respect to a by using the last expression of Eq. (JTTJ) to 
obtain 

To = - f r 2 - A r 3, (A15) 

where Eqs. ()A12| and (fT2|) has been used. Therefore we have 

t 5N ™ 5 5 5 rita 

fNL =6^T = Tr-3-6 r - (M6) 

Similarly, in order to calculate gNL, we can substitute Eqs. (|A12j) and (|A15|) into Eq. ()A14j) to write 

N aa = (- t- ~ ¥■) N °> (A17) 



K (J 3(7 3(7 

and then make derivative with respect to a. Finally we can obtain 



54 N aaa 9 1 , . , s 
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